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Abstract 

Civil  defense  shelters  are  often  constructed  beneath  the  ground  to  provide  protection  against  blast  loadings.  Concrete  is  widely  used  as  the 
material  for  the  defense  layer  of  the  shelters.  This  paper  adopts  a  continuum  damage  model  of  brittle  media  to  numerically  investigate  the 
dynamic  fracture  and  attenuation  effect  of  perforated  concrete  defense  layer  on  stress  waves  from  planar  charge.  The  model  includes  damage 
accumulation  and  loading  rate-dependence,  and  has  been  succinctly  implemented  into  the  dynamic  finite  element  code,  LS-DYNA,  via  its  user 
defined  subroutine.  The  numerical  results  reveal  that  the  adopted  model  can  well  predict  the  tensile  damage  owing  to  the  reflection  of  pressure 
waves  from  cavities  and  free  boundaries.  Again,  the  elastoplastic  properties  of  concrete  play  significant  roles  in  the  stress-wave  attenuation  and 
the  peak  values  of  hydrostatic  pressure  beneath  a  circular  cavity  are  largely  reduced.  One  empirical  formula  is  finally  proposed  to  relate  the  decay 
factor  of  peak  hydrostatic  pressure  to  the  cavity  dimensions  and  relative  position. 

©  2006  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  diffraction  and  attenuation  of  stress- waves  are  important 
aspects  in  the  design  of  civil  defense  shelters  [18,4,22,16].  The 
development  of  modern  military  technology  greatly  enhances 
the  destructive  power  and  hit  rate  of  new  weapons,  thus 
posing  great  challenges  to  defense  engineering  nowadays.  Civil 
defense  shelters  are  widely  constructed  beneath  the  ground  to 
provide  protection  against  blast  loadings  [18,1].  Blast  waves 
due  to  planar  charge  detonation  evolve  into  the  propagation 
of  stress-waves  in  underground  media.  Such  waves  will  be 
diffracted  when  hitting  obstacles  such  as  cavity,  crevice  or 
other  media.  The  wave  stress  intensity  may  be  greatly  reduced 
beneath  these  obstacles.  This  reduction  of  stress  from  the  peak 
value  is  often  termed  as  attenuation  or  insulation  of  stress- 
waves  [16]. 

Civil  defense  shields  are  widely  used  to  reduce  or  even 
prevent  the  damage  and  destruction  of  underground  structures 
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from  blast  or  impact  loading  [17,9,11,22,3].  When  designing 
a  civil  defense  work,  the  attenuation  of  and  screening  against 
stress-waves  must  be  considered  and  new  stress-wave  migration 
systems  are  expected.  Some  experimental  tests  related  to 
this  problem  have  been  performed  by  the  authors,  and  some 
field  observations  have  been  conducted.  Fig.  1  is  the  full- 
scale  test  appearance,  and  Fig.  2  depicts  the  configuration 
of  the  experimental  defense  structure.  It  clearly  shows  that 
a  layered  defense  structure  typically  consists  of  three  layers 
[18,12]:  a  soil  cover  layer,  a  protection  layer  and  a  support 
layer.  Furthermore,  the  protection  layer  has  two  sub-layers:  a 
projectile  shelter  layer  and  a  stress  distribution  layer.  Generally 
speaking,  the  stress  distribution  layer  is  made  of  perforated 
concrete  or  buffer  materials  such  as  sand  and  geofoam.  Its 
major  function  is  to  redistribute  the  blast  loading  over  a  larger 
area.  However,  studies  on  the  function  of  the  stress  distribution 
layer  are  limited  so  far. 

Artificial  cavities  are  normally  incorporated  into  the  stress 
distribution  layer  to  dissipate  the  energy  of  stress  waves. 
The  stress-wave  attenuation  due  to  the  dynamic  fracture  and 
plastic  deformation  of  the  perforated  concrete  is  far  from  well 
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Fig.  1.  Appearance  of  the  field  test  for  defense  structure  under  planar  charge. 
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Fig.  2.  Schematic  of  the  defense  structure  in  the  test  (unit  in  cm). 

understood  [19,20].  What  is  more,  the  existence  of  cavities  in  a 
perforated  defense  layer  may  create  great  potential  of  tensile 
damage  above  cavities  due  to  the  lower  wave  impedance  of 
air  inside  the  cavities.  Besides,  the  potential  damages  due  to 
reflected  tensile  waves  near  free  boundaries  (non-transmissive) 
should  also  receive  attention.  So  far,  studies  on  such  phenomena 
are  rarely  found  in  literature. 

In  the  present  study,  a  continuum  damage  model  for 
brittle  media  [21,13]  is  adopted  and  successfully  implemented 
into  the  finite  element  software,  LS-DYNA,  through  its  user 
subroutine  capability  [14].  The  numerical  scheme  with  the 
model  parameters  are  then  calibrated  against  limited  field  test 
data.  The  tensile  damage  owing  to  the  cavity  and  free  boundary 
reflection  is  investigated  subsequently.  The  attenuation  effect 
of  a  single  circular  cavity  on  peak  hydrostatic  pressure  is  also 
explored  in  detail  based  on  the  damage  model.  Finally,  an 
empirical  formula  is  proposed  to  relate  the  decay  factor  for  peak 
hydrostatic  pressure  to  the  cavity  parameters. 

2.  Continuum  damage  model  for  concrete  medium 

Within  a  concrete  mass,  flaws  and  cracks  exist  even  prior 
to  any  loading.  Under  externally  tensile  stresses,  these  flaws  or 
cracks  will  grow  in  size  and  number,  causing  the  deterioration 


of  concrete  stiffness  and  strength.  In  this  section,  a  continuum 
damage  model  for  brittle  concrete  was  adopted  and  a  software 
subroutine  was  specially  developed  for  the  study  of  the  tensile 
damage  and  attenuation  effects  of  perforated  concrete  defense 
layer  on  stress  waves. 

To  make  a  complex  problem  more  tractable,  it  is  assumed 
that  the  concrete  is  an  isotropic,  homogeneous,  continuous 
and  brittle  material  with  pre-existing  micro  cracks.  The  finite 
element  program  coupled  with  the  defined  model  captures  the 
general  effects  of  the  microcrack  system  on  material  stiffness, 
instead  of  treating  individual  cracks  which  are  direction- 
sensitive.  As  for  the  damage  process,  it  is  assumed  that 
damage  is  caused  by  the  activation  and  growth  of  pre-existing 
microcracks.  It  is  accumulated  over  time  and  is  irreversible. 

According  to  the  concepts  of  continuum  mechanics,  for 
isotropic  materials,  when  a  material  point  is  subjected  to 
stresses,  it  changes  in  volume  due  to  the  volumetric  parts 
of  the  stresses  and  in  shape  due  to  its  deviatoric  parts.  The 
volumetric  strain  sv  is  regarded  as  a  time-dependent  variable, 
which  determines  whether  the  microcracks  will  be  activated  and 
will  evolve. 

Under  external  loadings,  the  dynamic  fracture  of  concrete 
does  not  occur  unless  the  stress  is  larger  than  its  tensile  strength. 
This  is  accounted  for  by  setting  a  critical  value  sc  for  the 
volumetric  strain  sv  [15,13]: 

l-2v 

Sc  =  C Ts  (1) 

E 

where  E  is  the  Young’s  modulus,  v  is  the  Poisson’s  ratio  and  crs 
represents  the  tensile  strength  of  concrete. 

Obviously,  microcracks  will  be  activated  when  ev  exceeds 
sc,  whereas  the  microcrack  system  remains  stable  when  sv  is 
less  than  ec.  On  the  other  hand,  the  concrete  material  will  not 
fail  if  the  duration  of  an  applied  stress  whose  value  is  well  above 
the  tensile  strength  is  too  short.  In  this  regard,  a  crack  density 
Cd  is  defined  as  [21,13]: 

Cd  =  a (sv  -  sc)Pt  (2) 

where  a  and  are  material  constants,  t  is  time. 

The  differential  form  for  the  crack  density  is  written  as: 

dQ  =  afi(s  -  sc7~xtdsv  +  a(sv  -  sc)^dt.  (3) 

In  the  present  model,  the  damage  scalar  D  is  assumed  to  take 
the  following  form: 

D=  l-exp(-Cj).  (4) 

The  degraded  Young’s  modulus  Ed  for  the  damaged 
concrete  is  computed  as: 

Ed  =  E  •  (1  -  D).  (5) 

According  to  Jaeger  and  Cook  [8],  if  a  brittle  body  (rock, 
concrete)  contains  closed  cracks,  the  effective  Poisson’s  ratio 
is  greater  than  its  intrinsic  Poisson’s  ratio.  However,  when 
the  body  contains  very  flat,  open  cracks  or  equidimensional 
cavities,  the  effective  Poisson’s  ratio  is  smaller  than  its  intrinsic 
Poisson’s  ratio.  By  assuming  that  all  types  of  cracks  are 
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generated  equally  during  blast  damage,  the  Poisson’s  ratio  can 
be  assumed  to  remain  unchanged  for  the  damaged  body.  The 
use  of  an  unchanged  Poisson’s  ratio  also  has  the  advantage  that 
relationships  among  elastic  parameters  would  not  change  for 
the  damaged  body  [21]. 

As  for  a  concrete  mass,  when  it  is  under  tension  and  sv 
exceeds  sc,  the  incremental  relations  for  stress-strain  is: 


A  Oij  —  Kd8ij  +  2Gd  Aeij  (6) 

where  Kj  and  G  j  are,  respectively,  the  degraded  bulk  and  shear 
moduli,  8ij  is  the  Kronecker  delta,  A oyy  is  stress  increment, 
A Skk  and  A eij  are  increments  of  volumetric  and  deviatoric 
strains,  respectively. 

When  total  volumetric  strains  are  less  than  the  critical 
volumetric  strain,  elasto-plastic  model  is  used  to  model  the 
stress-strain  properties  of  concrete  and  the  Von  Mises  yield 
criterion  with  kinematic  hardening  is  used  for  the  plasticity 
model  [15].  In  this  situation,  elastic  parameters  of  the 
undamaged  concrete  are  adopted  since  it  is  assumed  that  micro 
cracks  remain  closed. 

The  yield  condition  [6]  is  expressed  as: 

<P  =  A2  -a2  (7) 

where  A  =  f\Cj  —  a >ij)(sij  —  a>ij),  cry  denotes  the  current 
radius  of  yield  surface,  and  deviatoric  stress  Sij  is  calculated  as: 

1 

sij  —  &ij  —  “  &kk $i j  •  (8) 

The  center  of  the  yield  surface  in  deviatoric  stress  space  is 
denoted  by  a  moving  tensor  co ,  whose  initial  value  is  zero,  as 
follows: 

Acoij  =  ^(l-ir)Epe[jAt  (9) 

where  ifr  is  the  hardening  parameter  (zero  for  the  kinematic 
case),  At  denotes  the  time  step,  Ep  is  the  plastic  hardening 
modulus,  the  current  radius  of  yield  surface  ay  is  the  sum  of 
initial  yield  strength  ao  and  the  growth  xf/EpS^: 

cry  =a0  +  fEpels  (10) 

where  the  effective  plastic  strain  is  defined  by: 

The  plastic  strain  ratio  is  the  difference  between  total  and 
elastic  strain  (right  superscript  e)  ratios: 


£P. 

ij 


:  Su  ~  £ 


ij’ 


(12) 


In  Eq.  (7),  (p  <  0  is  for  elastic  or  neutral  loading  while 
<p  >  0  is  for  plastic  hardening.  The  effective  plastic  strain  £P{{ 
and  stress  tensors  Oij  at  time  tn+\  respectively  are  as  follows: 


(£eff)n+1  —  (^) 

=  s"+1  +  pn+l8ij  (14) 


where, 


V/+1  =  Sfj  -  76GA<ff 
eff  3  G  +  Ep 


(15) 

(16) 


in  which,  <r"  is  the  yield  limit  at  time  tn,  S*j  and  pn+l  represent 
the  elastic  deviatoric  stress  state  and  the  hydrostatic  pressure  at 
time  tn+ 1. 


3.  Establishment  and  calibration  of  numerical  scheme 


3.1.  Numerical  tool 


In  the  present  study,  the  authors  have  successfully 
implemented  the  continuum  damage  concrete  model  described 
above  into  the  software  package,  LS-DYNA,  as  a  constitutive 
augment.  This  LS-DYNA  software  from  the  Livermore 
Software  Technology  Corporation  can  handle  either  two- 
or  three-dimensional  problems  for  high-speed  impact  and 
explosion  [6,14]. 

3.2.  Determination  of  model  parameters 

As  is  known,  concrete  C30  is  usually  used  to  construct 
the  protection  layer  in  civil  defense  engineering.  The  typical 
physical  parameters  are  as  follows:  initial  density  po  = 
2240.0  kg/m3,  the  Young’s  modulus  E  =  30.0  GPa,  the 
Poisson’s  ration  v  =  0.18,  uniaxial  compression  strength  oc  = 
30.0  MPa.  Assuming  the  tensile  strength  to  be  one-tenth  of 
the  uniaxial  compression  strength,  the  critical  volumetric  strain 
sc  (=6.4  x  10-5)  can  be  calculated  using  Eq.  (1). 

In  the  continuum  damage  constitutive  model,  a  critical 
parameter  is  the  crack  density  as  defined  by  Eq.  (2),  which 
in  turn  is  greatly  influenced  by  the  parameters  a  and  f. 
Thereinto,  f3  is  the  material  constant  which  determines  the  rate- 
dependence  of  fracture  stress.  It  has  been  observed  by  many 
researchers  that  the  fracture  stress  for  most  brittle  materials  like 
rock  and  concrete  depends  on  the  cube  root  of  the  strain  rate 
and  can  be  taken  as  2.0  [13,21].  So  far,  not  much  experience 
is  available  on  the  evaluation  of  the  parameter  a  which  should 
normally  be  determined  by  test  data  (see,  for  example,  [13]). 
In  view  of  this  and  with  the  aim  to  validate  the  developed 
numerical  subroutine  for  the  continuum  damage  constitutive 
model,  a  field  blast  test  is  conducted  in  the  present  study  [12].  In 
the  test,  planar  charge  is  adopted  as  shown  in  Fig.  3  which  gives 
the  arrangement  of  the  blasting  fuses  in  the  shallow  ground 
level.  Since  the  rise  time  of  blast  loading  is  typically  very 
short  and  the  overpressure  decays  exponentially  after  the  peak 
value,  an  equivalent  triangular  pressure  pulse  (see  Fig.  4)  is 
used  to  approximate  the  blast  loading  which  is  applied  on  the 
projectile  shelter  layer  [2,5,20].  The  intensity  and  duration  of 
this  overpressure  depend  on  the  explosive  charge  weight  and 
the  stand-off  distance  which  usually  follow  the  cubic  laws  [10, 
2] .  In  the  present  study,  the  peak  value  Qm ax  is  taken  as  about 
120.0  MPa  and  occurs  at  time  tmSLX  =  0.0  ps  (micro  second). 
This  peak  overpressure  then  decreases  linearly  to  zero  at  time 
tc  =  100.0  ps. 
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Fig.  3.  Field  assembly  of  blasting  fuses  for  planar  charge. 


Q 


Fig.  4.  Triangular  impulse  loading  for  explosion. 

The  subsurface  configuration  has  been  shown  in  Fig.  2 
earlier,  which  also  serves  as  the  domain  of  finite  element 
simulations.  The  FE  simulations  feature  4020  quadrilateral 
elements,  inclusive  of  the  steel  plate  (see  Fig.  2)  which  is 
divided  into  1 80  elements  with  the  properties  listed  in  Table  1 .  It 
is  calibrated  that,  when  a  is  taken  to  be  8.5  x  1010,  the  numerical 
results  are  in  fairly  good  agreement  with  the  field  recorded  data. 
For  example,  Fig.  5  is  the  history  of  hydrostatic  pressure  at 
the  stress  gauge  A 2  (refer  to  Fig.  2)  which  is  on  the  symmetry 
axis  and  159.0  cm  below  the  lower  surface  of  the  steel  plate. 
It  is  clear  that  the  numerical  simulations  using  the  FE  code 


Fig.  5.  Comparison  of  hydrostatic  pressure  between  the  numerical  calculations 
and  the  field  data  from  stress  gauge  A^. 


Table  1 

Material  parameters  for  steel  plate  used 


P0  (kg/m3) 

E  (GPa) 

V 

Gy  (MPa) 

Et  (GPa) 

7900.0 

207.0 

0.29 

274.0 

2.1 

LS-DYNA  with  the  continuum  damage  model  match  favorably 
with  the  field  measured  data.  The  FEM  simulation  also  reveals 
tensile  damage  at  t  —  1000.0  [is  near  the  roof  above  the  cavities 
as  shown  in  Fig.  6.  Such  numerical  revelation  is  corroborated  by 
the  post-blast  field  observation  of  the  multiple  tensile  cracks  on 
the  roof  layer  as  shown  in  Fig.  7.  It  should  be  noted  that  due  to 
the  very  limited  field  data  available,  more  extensive  comparison 
of  numerical  predictions  and  field  test  data  is  not  available. 
Details  of  the  field  blast  test  can  be  referred  to  in  [12]. 

4.  Numerical  investigations 

4.1.  Definition  of  decay  factor 

The  established  numerical  scheme  encompassing  the 
developed  user  subroutine  of  the  continuum  damage  model 
was  used  for  further  investigations  on  the  attenuation  effects 


LS-DYNA  user  input 
Time  =  1000 


Fringe  Levels 
l.OOOc+OO-l 
9.000e-0 1  - 

8.000C-01  -I 

7.000e-01  - 
6.000e-01  -I 
5.000e-01  -I 
4.000e-01  - 
3.000e-01  - 
2.000e-01  - 
1.000e-01  - 
O.OOOc+OO  - 


Fig.  6.  Contours  of  tensile  damage  in  problem  domain. 
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Fig.  7.  Post-blast  tensile  damage  of  a  roof  of  one  projectile  shelter  layer. 

of  concrete  defense  layer  on  blast  stress-waves  in  the  case  of 
single  cavity.  To  evaluate  the  attenuation  extent  of  stress-waves 
due  to  a  cavity  in  an  intelligent  defense  layer,  a  dimensionless 
variable,  ‘decay  factor’  DF ,  is  defined  as  follows: 


where  ao  denotes  the  peak  value  of  stress  component  such  as 
crx,  <jy,  or  P  at  a  specified  position  without  cavity  in  the  defense 
layer,  and  a  is  the  peak  value  of  the  same  component  at  the 
same  position  with  a  cavity  in  the  defense  layer.  Eq.  (17)  reveals 
that  a  larger  DF  would  imply  a  superior  attenuation  effect. 

4.2.  Computational  domain  with  cavity 

Fig.  8  shows  the  simulation  of  stress  distribution  layer 
embedded  with  a  single  circular  cavity.  The  computational 
domain  for  this  two-dimensional  model  is  taken  as  600  cm  x 
300  cm  (width  x  height).  The  upper  apex  of  the  cavity  is  at 
a  depth  e  beneath  the  surface  and  the  circular  cavity  has  a 
radius  R.  Thus  the  coordinates  of  the  center  is  (0,  e  +  R).  Other 
dimensions  are  labeled  in  Fig.  8.  Introducing  the  term  relative 
distance  y  as  the  distance  between  point  A  and  the  lower  apex  of 
the  cavity,  the  dimension  Y  follows  the  following  relationship: 

Y  =  y  +  2R  +  e.  (18) 

Only  the  right  half  of  the  domain  is  selected  in  the 
computations  owing  to  symmetry.  As  shown  in  Fig.  8, 
the  right  domain  has  the  following  boundary  conditions: 
pressure  boundary  on  the  top,  non-transmission  boundaries 
(free  reflection)  on  the  right  and  lower  sides,  and  symmetric 
boundary  on  the  left  side. 

4.3.  Effect  of  cavity  size 

Five  radii  are  chosen  in  the  computations:  0.0,  0.15,  0.30, 
0.45  and  0.60  m,  respectively.  The  parameter  e  is  fixed  as 
0.30  m.  The  model  parameters  for  concrete  are  the  same  as 
stated  in  Section  3.  For  comparison  purpose,  computations 
are  also  carried  out  using  the  primitive  linear  elastic  (F-E) 
model  of  concrete.  The  variations  of  peak  hydrostatic  pressure 
P  with  relative  distance  y  (refer  to  definition  in  Fig.  8)  are 


Fig.  9.  Effect  of  different  radii  R  of  circular  cavity  on  peak  hydrostatic  pressure 
for  (a)  L-E  model;  and  (b)  continuum  damage  model. 

shown  in  Fig.  9(a)  for  the  F-E  model  and  Fig.  9(b)  for  the 
continuum  damage  model.  In  both  cases,  the  peak  P  decreases 
with  increasing  radius  R  for  a  specified  position  y.  These  results 
reveal  that  a  larger  cavity  is  more  efficient  in  dissipating  the 
wave  energy.  On  the  other  hand,  the  peak  P  increases  with 
relative  distance  y  for  a  nonzero  R ,  which  indicates  that  wave 
energy  experiences  the  maximum  attenuation  around  the  lower 
apex  of  the  cavity  and  less  attenuation  further  down  below  the 
cavity.  Comparing  the  magnitudes  of  the  peak  P  in  Fig.  9(a)  and 
(b),  it  can  be  seen  that  the  elastoplastic  properties  of  concrete 
plays  a  crucial  role  on  the  magnitude  of  peak  P .  For  instance, 
when  R  =  0.30  m  and  y  =  0.60  m,  the  peak  P  is  30.4  MPa  for 
the  F-E  model  and  only  13.6  MPa  for  the  damage  model.  The 
change  of  decay  factor  DF  with  R  is  shown  in  Fig.  10  for  the 
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Fig.  10.  Variation  of  decay  factors  for  peak  P  with  different  radii  R  of  cavity 
for  continuum  damage  model. 

damage  model.  It  clearly  shows  that  DF  increases  with  R  but 
reduces  with  increasing  relative  distance  y. 

Fig.  9  also  shows  that  the  curves  for  the  case  of  zero  R 
(i.e.  without  cavity)  are  located  at  the  highest  positions  of  the 
plots,  which  demonstrates  the  effectiveness  of  a  cavity  in  the 
attenuation  of  stress  waves.  The  hydrostatic  pressure  contours 
for  cases  with  and  without  cavity  are  further  demonstrated 


in  Fig.  11.  Fig.  11(a)  shows  that  the  stress-waves  are  greatly 
disrupted  due  to  the  presence  of  the  cavity.  A  tensile  stress 
region  is  observed  above  the  cavity,  with  a  maximum  value  of 
almost  4.6  MPa.  For  a  concrete  grade  C30,  such  large  tensile 
stress  would  obviously  exceed  its  tensile  strength  and  lead  to 
the  tensile  damage  on  the  roof  above  the  cavity.  A  zero-stress 
region  can  be  noted  in  the  adjacent  domain  beneath  the  cavity. 
The  peak  P  beneath  the  cavity  is  evidently  lower.  In  addition, 
the  reflected  tensile  waves  can  also  be  observed  near  the  left  and 
right  boundaries.  For  the  non-cavity  case  shown  in  Fig.  11(b), 
the  tensile  stress  zone  is  only  located  at  the  side  boundaries 
of  the  calculation  domain  and  the  continuity  of  wave  front  is 
hardly  interrupted. 

Fig.  12  compares  the  damage  contours  for  a  cavity  radius 
R  =  0.6  m  at  t  =  440  [is  for  the  L-E  model  (Fig.  12(a))  and  the 
present  damage  model  (Fig.  12(c)).  As  an  added  comparison, 
the  computation  using  the  well-known  Johnson-Holmquist- 
Concrete  (J-H-C)  model  [7]  is  also  conducted  and  shown  in 
Fig.  12(b).  It  can  be  seen  that  the  L-E  model  is  not  capable 
of  predicting  damage  at  all.  The  J-H-C  model  is  capable  of 
evaluating  the  compression  damage  near  the  upper  surface  [20], 
but  is  incompetent  to  calculate  the  tensile  damage.  It  is  noted 
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Fig.  11.  Contours  of  hydrostatic  pressure  for  the  case  (a)  with  cavity  ( R  =  0.3  m  and  e  =  0.6  m);  and  (b)  without  cavity  (unit:  Mbar.  1  Mbar  =  lO*5  MPa). 
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Fig.  12.  Contours  of  damage  for  R  =  0.60  m  cavity  at  t  =  440  \is  under  the  same  blast  loading  for  (a)  L-E  model;  (b)  J-H-C  model;  and  (c)  continuum  damage 
model. 


that  the  proposed  damage  model  used  in  the  present  study  can 
well  evaluate  the  tensile  damage  above  the  cavity  and  near 
the  side  boundaries.  In  particular,  the  spalling  at  those  places 
may  occur  with  the  increase  of  reflected  wave’s  magnitude. 
Hence,  special  protective  measures  at  these  identified  locations 
are  highly  desirable  when  the  applied  blast  loading  is  very  large. 

4.4.  Effect  of  cavity  position 

The  effect  of  the  cavity  position  e  on  the  attenuation  of 
pressure  waves  is  shown  in  Fig.  13.  When  R  is  fixed,  the 
decay  factors  DF  are  practically  identical  for  all  e  values.  This 


phenomenon  can  also  be  observed  from  the  time-hydrostatic 
pressure  curves  at  y  =  0.5  m  shown  in  Fig.  14.  The  peak  value 
of  hydrostatic  pressure  does  not  vary  much  for  all  e  values. 
These  numerical  results  reveal  that  the  overhead  cover  of  the 
cavity  in  the  stress  distribution  layer  may  not  be  a  critical  factor 
in  the  design  considerations  of  a  concrete  defense  layer. 

5.  Empirical  formula  for  decay  factor 

Based  on  the  above  numerical  results,  it  is  established  that 
DF  is  essentially  a  function  of  R  and  y: 

DF  =  f(R,y) 


(19) 
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Fig.  13.  Variation  of  decay  factors  for  peak  hydrostatic  pressure  with  different 
positions  e  of  cavity. 


Fig.  15.  Comparison  of  empirical  formulae  and  numerical  results  of  decay 
factors  for  the  peak  hydrostatic  pressure. 


Fig.  14.  Effect  of  cavity  position  e  on  hydrostatic  pressure  history  at  y  =  0.5  m 
point  on  7 -axis. 


where  the  relative  distance  y  =  Y  —  2R  —  e. 

In  order  to  determine  the  particular  form  of  Eq.  (19),  one 
value  of  e  (0.30  m),  four  values  of  R  (0.15,  0.30,  0.45  and 
0.60  m)  and  five  values  of  y  (0.20,  0.40,  0.60,  0.80  and  1.0  m) 
are  chosen  from  Fig.  10.  Through  a  curve  fitting  regression 
analysis,  DF  for  the  peak  P  at  a  specified  point  on  the  7-axis  is 
approximately  expressed  as  follows: 


DF  =  exp 


(20) 


This  formula  shows  that  DF  increases  with  R  and  decreases 
with  y.  When  the  term  ^  approaches  zero,  DF  is  almost  equal 
to  unity.  Fig.  15  shows  the  comparison  of  numerical  results  with 
the  predictions  from  the  above  empirical  formula.  It  is  clear  that 
this  empirical  formula  agrees  favorably  well  with  the  numerical 
results.  Therefore,  Eq.  (20)  provides  a  very  useful  and  practical 
means  to  expeditiously  estimate  the  required  size  of  a  circular 
cavity  if  the  expected  DF  (at  a  relative  distance  y)  and  position 
e  are  known.  On  the  other  hand,  when  DF  and  R  are  given, 
the  relative  distance  y  can  also  be  determined  readily  from 
Eq.  (20). 


6.  Conclusions 


A  numerical  scheme  encompassing  a  developed  FEM 
subroutine  for  the  continuum  damage  concrete  model  was 


established  in  the  present  study  and  calibrated  against  field 
blast  tests.  The  established  numerical  tool  was  then  used  for 
further  investigation  of  the  tensile  damage  and  attenuation 
effect  of  concrete  defense  layer  on  blast  stress  waves.  From 
these  studies,  the  following  conclusions  and  understandings  can 
be  drawn. 

(1)  Compared  to  other  models  like  the  L-E  model  or  J-H-C 
model,  the  continuum  damage  concrete  model  presented 
in  this  study  can  not  only  capture  the  tensile  damage, 
but  can  also  evaluate  the  material  plastic  behaviors  during 
compression.  The  present  numerical  scheme  incorporating 
this  damage  model  possesses  great  potential  and  advantage 
for  future  analyses  and  design  of  perforated  concrete 
defense  layer  against  blast-induced  stress  waves. 

(2)  Artificial  cavities  embedded  in  the  concrete  defense  layer 
have  a  vital  effect  on  the  attenuation  of  stress- waves.  As  for 
a  circular  cavity,  the  radius  R  has  a  significant  effect  on  the 
reduction  of  stress-waves  according  to  our  computations, 
while  the  effect  of  the  position  of  the  cavity  is  found  to  be 
negligible. 

(3)  One  empirical  formula  is  obtained  to  correlate  the  decay 
factor,  dimensions  of  a  circular  cavity  and  relative  distance. 
The  formula  can  be  employed  to  determine  the  minimum 
radius  of  a  cavity  Rm in,  or  minimum  safe  distance  ymjn 
when  the  decay  factor  is  given.  Therefore,  it  can  serve  as 
a  very  practical  and  handy  tool  in  the  process  of  designing 
an  intelligent  defense  layer  in  civil  defense  engineering. 

It  should  be  noted  that  the  third  conclusion  is  applicable 
only  when  the  damage  of  concrete  does  not  initiate  the  collapse 
of  the  cavity.  Further  work  is  necessary  when  the  damage  of 
concrete  defense  layer  is  severe  enough  to  cause  a  complete 
collapse  of  the  cavity. 
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